Enhanced polarization of CMB from thermal gravitational waves 
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If inflation was preceded by a radiation era then at the time of inflation there will exist a decoupled 
thermal distribution of gravitons. Gravitational waves generated during inflation will be amplified by 
the process of stimulated emission into the existing thermal distribution of gravitons. Consequently 
the usual zero temperature scale invariant tensor spectrum is modified by a temperature dependent 
factor. This thermal correction factor amplify the B-mode polarization of the CMB by an order of 
magnitude at large angles, which may now be in the range of observability of WMAP. 
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Inflation [1], in addition to solving the horizon and 
flatness problems of the standard hot big-bang model, 
generates a nearly scale invariant density perturbations, 
which has been tested in the observations of the CMBR 
angular spectrum. One prediction of inflationary mod- 
els which has not yet been tested is the existence of a 
nearly scale invariant spectrum of gravitational waves 
0, [1| • The definitive test of the existence of these cos- 
mological gravitational waves would be the observation 
of B mode polarization in the CMB. The recent WMAP 
three year results [3| give only an upper bound on the B 
mode polarization, ^CS_ 6) < 0.05(,uif) 2 . 



In this paper we show that if inflation was preceded by 
a radiation era then there would be a thermal background 
of gravitons at the time of inflation. This thermal distri- 
bution of gravitons would have decoupled close to Plank 
era. The generation of tensor perturbation during infla- 
tion would be by stimulated emission into this existing 
thermal background of gravitational waves. This process 
changes the power spectrum of tensor modes by an extra 
temperature dependent factor coth(fc/2T). At large an- 
gular scales (I < 30) the power spectrum Pt — At k nT 
of gravitational waves generated during inflation would 
have a spectral index = — I — 2e, instead of the stan- 
dard slow roll inflation prediction n-r — — 2e which im- 
plies that C^f ~ 10 x Cj?.^. If a thermal enhancement 
of low / BB modes exists then it will be observable with 
WMAP or in the upcoming Planck [j| experiment. In 
the conclusion, we discuss the implications of inflation- 
ary models from the observations or non-observation of 
this low I thermal enhancement. 



The tensor perturbations have two independent de- 
grees of freedom which can be chosen as h + and h x polar- 
ization modes. To compute the spectrum of gravitational 
waves h(x, r) during inflation, we express and 



in terms of the creation- annihilation operator as 
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where a(r) is the scale factor, k is the comoving 
wavenumber, k = |k|, and M p = 1.22 x 10 19 GeV is the 
Plank mass and i = +, x. The power spectrum of the 
tensor perturbations is defined as: 

{h k h k ,) = ^P T 6 3 (k~k r ). (2) 

The usual quantization condition between the fields and 
their canonical momenta yields [a k , aL] = <5 3 (k — k') and 
the vacuum satisfies ak|0) = 0. If the graviton field 
had zero occupation prior to inflation then [a k , a fc ,] = 
<5 3 (k — k') and the vacuum satisfies ctk|0) = 0. If the 
graviton field had zero occupation prior to inflation then 
(a k a k ) = and we would obtain a correlation function 
~ |//c(t)| 2 . However if the graviton field was in thermal 
equilibrium at some earlier epoch it will retain its thermal 
distribution even after decoupling from the other radia- 
tion fields and its occupation number will be given by: 
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Using Eq. ([1]) and Eq. ((3|) it can be seen that 
16^|/ fc (r)| 2 



(3) 



a 2 (r)M 2 

16^|/ fc (r)| 2 
a 2 (r)M 2 



1 



coth 



2 


k 


- 1 




" k ' 


5 


2T 





5 3 (k-k') 



j (k-k') 



(4) 



From the defining relation, Eq. ([2]), for the tensor power 
spectrum and Eq. (|4|) we find that the power spectrum 
for the thermal inflatons can be expressed in terms of the 
mode functions /fe(r) as: 
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The mode functions /fe(r) obey the minimally coupled 
Klein-Gordon equation: 



a 
a 
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In a quasi De Sitter universe during inflation, conformal 
time r (dr = dt/a) and the scale factor during inflation 
a(r) are related by a(r) = —1/Hr{l — e) where e = 

Yg^r I J and V is the potential for the inflaton field. 

For constant e the mode functions fk(r) obey the min- 
imally coupled Klein-Gordon equation 



fk = 0, 



where k = |k| and, for small e and <5, v 
has the general solution given by, 
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When the modes are well within the horizon they can 
be approximated by flat space-time solutions /fe (t) = 
-j=e~ lkT , (fc ^ ai/). Matching the general solution in 
Eq. ([5]) with the solution in the high frequency (flat space- 
time) limit gives the value of the constants of integration 
Cl (k) = ^e'^+^f and c 2 {k) = 0. Eq. © then 
implies that for — kr > 1 or fc < aH, 
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Substituting the solution as given in Eq. ([9]) for the super- 
horizon modes (k <C aH) in Eq. ([5]) for the tensor power 
spectrum, we obtain: 
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with tit = 3 — 2v = —2e. We can now rewrite the power 
spectrum as, 



P T (k) = A T (k ) 
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where k is referred to as the pivot point and A(kg) is 

the normalization constant, At^q) — jfe (y-^j-^j where 
Hk is the Hubble parameter evaluated when aH = ka 
during inflation. 

The angular power spectrum of the BB polarization 
modes generated by the gravitational waves is given by 
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where g{rf) — he~ K is the visibility function and k is the 
differential optical depth for Thomson scattering. The 
EE polarization signal gets a contribution from the ten- 
sor perturbations but it is mainly dominated by the scalar 
perturbations. So the best signal for gravitational waves 
is the BB polarization angular spectrum which is gener- 
ated by the primordial tensor perturbations only. 

The temperature dependent factor becomes important 
when the ratio k/(2T) is less than unity. The co-moving 
wave-number k and the co-moving temperature T can be 
related to the physical parameters at the time of inflation 
as follows. Taking the largest measurable perturbation 
scale k now /a now ~ R^ 1 (where Rh — 4000Mpc is the size 
of the present horizon) , and assuming that perturbations 
of the present horizon scale were just leaving the infla- 
tionary horizon H^ 1 at the beginning of inflation we see 
that the temperature at the beginning of inflation Tij ai 
must be, 
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in order to have a significant effect on the tensor power 
spectrum. Ti/ai ~ (30V/ g^n 2 ) 1 / 4 , V being the inflaton 
potential which is related to the curvature at the time of 
inflation, H = (87r/3) 1 / 2 y 1 / 2 /M P and .g* - 100 is the 
effective number of spin/polarization degrees of freedom 
of relativistic particles. Therefore inflation is expected to 
start at a temperature Ti/a t = 0.24(iJMp) 1 / 2 . Actually 
the gravitons which are decoupled will have a temper- 
ature slightly below the radiation temperature because 
of the particles (like the inflaton itself) which have an- 
nihilated into radiation prior to inflation. But as the 
effective number of degrees of freedom, ~ 100, is large 
this difference of temperature is not significant. So for in- 
flation at the GUT scale, V 1 / 4 - 10 15 GeV, we have H ~ 
10 11 GeV and the temperature at the start of inflation 
Ti/ai ~ 10 14 ~ 15 GeV. So the enhancement of the gravi- 
ton power spectrum by the factor coth^^) = coth(-^ i ) 
could be by as large as a factor of 10 4-5 at low k due to 
thermal gravitons. 

In Fig. [1] we show the angular correlations of CMBR 
temperature and polarization assuming a thermal gravi- 
ton spectrum ( along with the WMAP three years data 
0]). The plots for TT, TE and BB correspond to co- 
moving graviton temperature T = .OOlMpc -1 . For com- 
parison we have plotted the BB angular correlations at 
T = 0. We see that with a temperature T = .OOlMpc -1 
the BB correlations are amplified at I < 30. We see that 
only the BB correlation is enhanced by the correction to 
the tensor power spectrum as expected. The contribution 
of tensors to the TT angular spectrum is comparable at 
low I to the contribution from the scalars and there exists 
the possibility that this large tensor contribution at low I 
may be detected from the analysis of the the TT angular 
spectrum alone. 
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We have added the unlensed scalar and tensor con- 
tributions to generate the TT, EE, TE and BB correla- 
tions. The plots were obtained by running CMBFAST 
with the following parameters fif, = 0.05 , f2 c = 0.25 
and Q v = 0.70. Optical depth r = 0.08 and Hubble 
parameter h = 0.7. The value of scalar spectral index 
n s = 0.97 and the value of tensor spectral index is taken 
tit = —0.01. Tensor to scalar ratio is taken to be r(ko) = 
0.1 at fc = 0.002Mpc~ x . The output of the CMBFAST 
was normalized to the WMAP values at k — 0.002Mpc _1 
(i.e / = 30). For the curves shown in Fig. [T] the tensor 
power spectra is modified due to thermal effects with 
^ = 500fc. At k = 0.0002MPC" 1 , coth(500fc) = 10 so 
there is a large enhancement of the BB polarization at 
I = 2 - 6, while at k = 0.002Mpc _1 , coth(500fc ) ~ 1.3 
and there is hardly any enhancement of the BB signal (or 
in the value of r(ko) in keeping with the observational 
constraints from WMAP+SDSS [Icj]). The magnitude of 
the co-moving graviton temperature needed to produce 
this effect is T/a n0 w = 10 _3 Mpc _1 . This corresponds to 



a temperature of T /a,i ~ 4 x R h 



x a r , 



•j/di = 4iJ (where 



Rh ~ 4000Mpc is the size of the present horizon). As we 
have seen inflation can start as soon as the temperature 
Ti/cii falls below V 1 ! 4 ~ 10 4 _ff. Consequently a temper- 
ature larger than AH at the beginning of inflation is not 
unreasonably high. 

In standard inflation models the vacuum fluctuations 
of the inflaton field give the density perturbations. The 
inflaton is assumed to be decoupled from the radiation 
at the onset of inflation, however if there was a radia- 
tion era prior to inflation then the scalar curvature power 
spectrum is modified by the same temperature dependent 
factor [llj as the tensor power spectrum in Eq. (jTUJ) , 



Pn(k) = 



H 4 



4tt 2 



/ k N 


n»— 1 


' k ' 




| coth 






2T 



(14) 



The extra temperature dependent term implies that there 
should be an up-turn of the TT anisotropy spectrum 
at low I . This expected up-turn in 1(1 + 1 )(7; is not 
seen in the WMAP one-year TT spectrum [ll|. This 
means that there is no significant number density of back- 
ground density of inflatons at the time when the modes 
which are currently entering our horizon, were exiting 
the horizon during inflation. This could happen for two 
reasons. The background density of inflatons may have 
decayed or annihilated into lighter particles by this time 
or the inflaton was cooled from the expected tempera- 
ture of 0.24(ifMp) 1 / 2 to below H by the time the modes 
corresponding to our present horizon were leaving the 
De-Sitter horizon. This implies that there were an ex- 
tra AiV = ln(0.24(M P /i/) 1 / 2 ) e-foldings (which has the 
value AiV ~ 10 for GUT scale inflation) than what is 
needed to solve the horizon problem. In the case of gravi- 
tons the first condition does not apply as they decouple 
at the Planck scale and if the expected upturn in the 



BB mode spectrum is not seen that would imply that 
the duration of inflation was longer than what is needed 
to solve the horizon problem. 

In warm inflation models [13] where the inflaton is 
in thermal equilibrium with the radiation bath and the 
scalar curvature perturbations are generated by thermal 
fluctuations instead of by quantum fluctuations, there is 
no coth(fc/2T) correction in the inflaton power spectrum 
due to stimulated emission. However this correction fac- 
tor will be present in the graviton spectrum since gravi- 
tons are still produced by quantum fluctuations, in warm 
inflationary models. The scalar curvature perturbation 
in warm inflation is 11311: 
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where T designates the decay width of the inflaton field 
and T r is the temperature of the radiation bath. 
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FIG. 1: The TT, TE and the BB correlations with thermal 
graviton spectrum along with the WMAP three years data [4]. 
The plots for TT , TE and BB correspond to co-moving graviton 
temperature T = .OOlMpc -1 . For comparison we have plotted 
the BB angular correlations at T = 0. We see that with a 
graviton temperature T = .OOlMpc" 1 the BB correlations are 
amplified at I < 30. 

There are observational constraints on the tensor 
scalar ratio defined as: 
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From the combination of WMAP three year data [lOj 
and SDSS large scale structure surveys [l4j we have the 
bound r(k = O^Mpc" 1 ) < 0.28(95%CL) where k = 
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0.002Mpc~ corresponds to I = rofco — 30 with the dis- 
tance to the decoupling surface r = 14, 400Mpc. SDSS 
measures galaxy distributions at red-shifts a ~ 0.1 and 
probes k in the range 0.016/iMpc _1 < k < O.llWVlpc -1 . 
From the expressions of P-r. in warm inflation, Eq. (|15[) . 
and Pt we see that the scalar-tensor ratio in warm in- 
flation models (assuming a nearly scale invariant tensor 
power spectrum) has a scale dependence at large angles 
given by: 
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We see that r(k) has a spectral index tit ~ — 1 for large 
scale perturbations. If we consider fe ~ 0.0002Mpc _1 
which corresponds to I ~ 3 then the value of r(fc) = 
lOr(fco). So even with r(k ) ~ 0.1 as constrained by 
galaxy surveys, we can have r(k) ~ 1 at the quadrupole 
anisotropy. The -B mode polarization at Z = 3 is enhanced 
from its value at I = 30 by a corresponding factor of 10. 
This is true as long as the temperature Ti/ai < 10 4 H 
which as we have seen in the earlier discussion is expected 
if there is a thermal era prior to inflation. 

For example taking the inflaton potential to be V = 
(l/2)m 2 2 , we have the scalar power: 
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and the tensor power, 
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and the scalar-tensor ratio, 



r(fc ) = 8.413 
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-B-mode polarization in the CMB would confirm the ex- 
istence of primordial super-horizon gravitational waves. 
Observationally, the three year WMAP data only gives 
an upper bound on Cf B with I = (2 — 6) 0] . The error 
bars on the C BB are presently a factor of five larger than 
the predictions from standard inflation theory with scalar 
tensor ratio as large as 0.3, which is close to the observa- 
tional upper bound ro.002 < 0.28(95%CL). In this paper 
we show that due to thermal gravitons, the C BB at low 
I ~ (2 — 6) could be larger by a factor of 10 compared 
to what would be expected from the observational con- 
straint on r and could be within the range of observabil- 
ity of WMAP. The upcoming Planck experiment [j| will 
measure C^ 1 _ w -. at the level of 10~ 4 (fiK) 2 . Ground 
based polarization experiments [3] like QUaD, QUIET, 
Clover and PolarBear measure anisotropies at small an- 
gular scales only (at I > 100 where thermal effects dis- 
cussed in this paper are negligible) and can observe C BB 
at the level 10~ 2 (^K) 2 . These experiments can probe 
r in the range 0.05 — 0.1 independent of thermal effects. 
A combination of data from WMAP/Plank at large an- 
gles and ground based polarization experiments at small 
angles will therefore either observe or definitely rule out 
the thermal enhancement effect. 

If WMAP or Planck rule out a spectral index of 
ji t ~ -i at low I, which is the prediction from ther- 
mal gravitons, then for the standard inflationary mod- 
els it would mean that the duration of inflation has to 
be longer by AiV = ln(0.24(M P / J ff) 1 / 2 ) e-foldings than 
what is needed to solve the horizon problem. Warm in- 
flation models [13, US] cannot evade this constraints by 
supercooling during inflation. If _B-modes are observed 
and the tensor spectral index at low I is not close to — 1, 
then warm inflation models can be ruled out. 



where 4>o is the value of the inflaton field when the scale 
ko = 0.002Mpc _1 was leaving the inflaton horizon. By 
choosing the parameters m = 1.4 x 10 12 GeV, T — 0.5 x 
10 13 GeV, T ~ T r = 0.24 x 10 16 GeV, 4> Q ~ 0.8 x 10 19 GeV 
we have Pr ~ 2.3 x 10~ 9 as required by WMAP three 
year data and r(ko) = 0.095. The value of r is larger at 
k = 0.0002Mpc _1 by a factor of - 10 and the B-modes 
are magnified at / = 3 compared to their value at I = 30 
by a factor 10, also in warm inflation scenarios. 

Direct observation of gravitational waves would nail 
the last still unconfirmed prediction of inflation. The am- 
plitude of gravitational waves gives the Hubble curvature 
during inflation and would tell us the value of the infla- 
tion potential 15j]. In addition gravitational waves pro- 



duced during inflation can have several applications like 
leptogenesis by the gravitational spin-coupling to neu- 
trinos [HI ] or by a gravitational Chern-Simon coupling 
of the lepton number current [17] . Observation of the 
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